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We study the percolation phase transition in hierarchical scale-free nets. Depending on the method 
of construction, the nets can be fractal or small-world (the diameter grows either algebraically or log- 
arithmically with the net size), assortative or disassortative (a measure of the tendency of like-degree 
nodes to be connected to one another), or possess various degrees of clustering. The percolation 
phase transition can be analyzed exactly in all these cases, due to the self-similar structure of the hi- 
erarchical nets. We find different types of criticality, illustrating the crucial effect of other structural 
properties besides the scale-free degree distribution of the nets. 
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I. INTRODUCTION 

Many large complex nets, such as the Internet and the 
World Wide Web, social networks of contact, and net- 
works of interactions between proteins are scale- free: the 
degree ki (number of links attached to a node, i) has 
a distribution with a heavy power-law tail, P(k) ~ fc~ 7 . 
Because of their ubiquitousness in everyday life the struc- 
ture and physical properties of scale-free nets have at- 
tracted much recent attention 

The percolation problem is of particular practical in- 
terest: Is the integrity of the Internet compromised fol- 
lowing random breakdown of a fraction of its routers? 
What fraction of a population ought to be vaccinated to 
arrest the spread of an epidemics that spreads by social 
contact? Initial studies of percolation addressed the case 
of stochastic scale-free nets, where the links between the 
nodes are drawn at random, so as to satisfy the scale- 
free degree distribution (for example, by the algorithm 
due to Molloy and Reed [2|). These studies showed that 
scale-free nets are resilient to random dilution, provided 
that the degree exponent 7 is smaller than 3. Explicit 
expressions for the critical exponents characterizing the 
transition as a function of 7 were also derived [!, 0> S @] • 

Stochastic Molloy-Reed scale-free nets are limited, 
though. Having fixed the degree distribution, all other 
structural properties (such as the extent of clustering, as- 
sortativity, etc.) are fixed as well, in contrast with man- 
made and natural scale-free nets that show a great deal 
of variation in these other properties. In this paper, we 
study percolation in hierarchical scale- free nets @, H, Q ■ 
Hierarchical scale-free nets may be constructed that are 
small- world or not, and with various degrees of assorta- 
tivity, clustering, and other properties [10 ]. 

Hierarchical nets have been studied before, as exotic 
examples where renormalization group techniques yield 
exact results 0, H, Q, including the percolation phase 
transition and the q — * 1 limit of the Potts model [lll[l2|]. 



We study percolation directly, by focusing on the size of 
the giant component (the largest component left after 
dilution) and the probability of contact between hubs 
(nodes of highest degree). Our aim is to elucidate the 
effect of the various structural properties of the nets on 
the percolation phase transition. As we shall see below, 
whether the transition takes place or not, and its char- 
acter, depends not only on the degree exponent (as in 
stochastic nets) but also on other factors [131 ]. 



II. HIERARCHICAL SCALE-FREE NETS 

Hierarchical scale-free nets are constructed in a recur- 
sive fashion. We focus on the special class of (u, v)- 
flowers [Tf3 [. where each link in generation n is replaced 
by two parallel paths consisting of u and v links, to yield 
generation n+1. A natural choice for the genus at genera- 
tion n = 1 is a cycle graph (a ring) consisting of u+v = w 
links and nodes (other choices are possible). The case of 
u = 1, v = 2 (Fig. [T]) has been studied previously by 
Dorogotsev, Goltsev and Mendes (DGM) Q. In the 
following we assume that u < v, without loss of general- 
ity. 

All (u, v)-flowers are self-similar, as evident from an 
equivalent method of construction: to produce genera- 
tion n + 1, make w copies of the net in generation n and 
join them at the hubs (the nodes of highest degree), as 
illustrated in Fig. [TJ;. 

It is easy to see, from the second method of construc- 
tion, that the number of links (the size) of a (u, v)-flowcr 
of generation n is 



M n = (u + v) n = w n . 



(1) 



At the same time, the number of nodes (the order) obeys 
the recursion relation 

N n = wNn^x - w , 

which, together with the boundary condition N\ = w, 
yields 
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FIG. 1: The (1, 2)-flower, or DGM network, (a) Method of 
construction: each link in generation n is replaced by two 
paths of 1 and 2 links long, (b) Generations n = 1,2,3. 
(c) Alternative method of construction: generation n + 1 is 
obtained by joining three replicas of generation n at the hubs 
(marked by A,B,C). 



Similar considerations let us reproduce the full degree 
distribution. By construction, (it, v)-flowers have only 
nodes of degree k — 2 m , m = 1, 2, . . . , n. Let N n (m) be 
the number of nodes of degree 2 m in the (tt, v)-flower of 
generation n, then 

N n (m) = JVn_i(m -l) + (w- 2) w "- 1 <5 m ,i , 



leading to 



N n {m) 



(w - 2)w r ' 



, m < n, 
m = n. 



(3) 



As in the DGM case, this corresponds to a scale-free de- 
gree distribution, P(k) ~ fc -7 , of degree exponent 



7 = 1 



ln(it + v) 
ln2 



(4) 



The self-similarity of (it, u)-nets, coupled with the fact 
that different replicas meet at a single node, makes them 
amenable to exact analysis by renormalization group 
techniques. 



Network diameter and dimensionality 

There is a vast difference between (it, i>)-flowers with 
it = 1 and it > 1. If it = 1 the diameter L n of the n-th 
generation flower (the longest shortest path between any 
two nodes) scales linearly with n. For example, L n = n 
for the (1, 2)-flower Q and L n = 2n for the (1, 3)-fiower. 
It is easy to see that the diameter of the (1, v)-flower, for 
v odd, is L n = (v — l)n + (3 — v) /2, and, while deriving a 



similar result for v even is far from trivial, one can show 
that L n ~ (u — l)n. 

For it > 1, however, the diameter grows as a power of 
n. For example, for the (2, 2)-flower we find L n = 2™, 
and, more generally, if u + v is even (and u > 1), 

L n = — + hi" ] 



u- 1 



u- 1 



For it + v odd one may establish bounds showing that 
L n ~ it™. To summarize, 



(v — l)n it = 1, 
u n it > 1. 



(5) 



Since AT,, ~ (it + v) n , we can recast these relations as 
J lnA^ u = 1, , . 

^ ~ S jyInu/ln(u+«) u > 1. ( - D ' ) 

Thus, for it = 1 the flowers are small world, similar to 
stochastic scale-free nets with 7 > 3. For u > 1 the nets 
are in fact fractal, with fractal dimension 



d = 



ln(u + v) 
ln u 



it > 1 , 



(7) 



since the mass increases by it + v (from one generation 
to the next) while the diameter increases by u. (1, v)- 
flowers are m/tmie-dimensional. In [ 101 ] we showed how 
these nets may be characterized by a different measure 
of dimension that takes into account their small-world 
scaling. 

The difference between flowers with u — 1 and u > 1 is 
perhaps best exemplified by the (1,3)- vs. the (2,2)-flower 
(Fig. [2]). The nets have identical degree distributions, 
node for node, with degree exponent 7 = 3 — similar to 
the famed Barabasi-Albert model [lj| — but the (1,3)- 
flower is small world (or infinite-dimensional), while the 
(2, 2)-flower is a fractal of dimension d = 2. 



Other structural properties 

Upon varying u and v the hierarchical flowers acquire 
different structural properties. Consider, for example, 
their assortativity — the extent to which nodes of similar 
degree connect with one another [l6j]. In the (1, 3) -flower, 
nodes of degree 2 m and 2 m+1 are only one link apart, and 
the assortativity index is 0; while in the (2, 2)-flower the 
same nodes are 2 m_1 links apart, and its assortativity 
index tends to —1/2 (as A^ — > 00), indicating a high 
degree of disassortativity, and more in line with naturally 
occurring scale- free nets [l7l . [l8| . More generally, the 
degree of assortativity, r, is r — ► (v — 3)/2u, for u = 1, 
andr-> -2/(u + v), for u > 1, (as A^ -> 00) QUE!- 

Another property of interest is clustering, a measure of 
the likelihood for neighbors of a node to be neighbors of 
one another [2fj| . (it, w)-flowers with u > 1 have zero clus- 
tering: the neighbors of a node are never neighbors of 
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FIG. 2: (it, i;)-nowers with u + u = 4 (7 = 3). (a) Small world: 
u — 1 and « = 3. (b) Fractal: u = 2 and w = 2. The graphs 
may also be iterated by joining four replicas of generation n 
at the appropriate hubs. 
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FIG. 3: Decorated (2, 2)-flower. Construction method (top): 
Each link is replaced by two parallel paths of u and v links 
long, and an additional link (broken line) that does not get it- 
erated, (bottom): The decorated (2, 2)-flower for generations 
n = 1,2,3. 

construction one can generate scale-free nets with differ- 
ing structural properties. By changing one property at a 
time one can then hope to understand their effect on var- 
ious physical phenomena, such as the percolation phase 
transition. One can also construct nets that mimic every- 
day life networks as closely as possible. The decorated 
(2, 2)-flower, with its degree exponent 7 = 3, disassor- 
tativity, and high degree of clustering, is a reasonable 
candidate for the latter. 



one another. The DGM net (u = 1, v = 2) has clustering 
coefficient c w 0.7998, and c gets smaller with increasing 
v (or degree exponent 7) , quite in line with the clustering 
coefficient of everyday life scale-free nets. 



Decorated flowers 

So far we have seen hierarchical nets that are either 
fractal and disassortative (u > 1), or small world and as- 
sortative (u = 1 and v > 3). It is also possible to obtain 
hierarchical nets that are small world and disassortative 
at the same time. One way to do this is by constructing 
a fractal (it, w)-flower (11 > 1) and adding a link between 
opposite hubs at the end of each iteration step: the addi- 
tional link does not get iterated 0,11, H- Fig.[3]illustrates 
this procedure for the case of the (2, 2)-flower. 

The added link has a negligible effect on the degree 
distribution: the degree exponent still approaches 7 = 
1 + ln(tt + v)/ In 2, as n —* 00. On the other hand, it has 
a dramatic effect on the diameter of the net, which now 
grows linearly in n (logarithmically in N), making it a 
small world network. The nets remain strongly disassor- 
tative: for example, for the (2, 2)-flower the assortativity 
changes from —1/2 (without) to —9/26 with the addition 
of the non-iterated links. Finally, the added links have a 
dramatic effect on the clustering of the nets, which grows 
from zero to about 0.82008, for the (2, 2)-fiower |[. 

The main point is that by manipulating the method of 



III. THE PERCOLATION PHASE TRANSITION 

We now turn to the study of percolation in hierarchical 
scale-free nets. The recursive nature of the (it, i>)-fiowers, 
coupled with their finite ramification, make it possible to 
obtain an exact solution by a real-space renormalization 
group analysis, including the finite-size behavior around 
the transition point. 

Our plan is as follows. We first study percolation in 
fractal hierarchical nets. Having finite dimensionality 
they resemble regular and fractal lattices, and the perco- 
lation phase transition is similar to what is found there 
as well. We then study percolation in the (1, v) -flowers, 
which are small world, as most everyday life complex net- 
works. Unlike everyday life nets, the (1, w)-flowers have 
no percolation phase transition, even for v > 3, or 7 > 3. 
Clearly, the (1, u)-fiowers fail to mimic everyday life net- 
works in some crucial aspect — perhaps their high as- 
sortativity. We therefore conclude with an analysis of 
the decorated (2, 2)-flower. The transition there most 
closely resembles that of everyday life nets, but some dif- 
ferences remain. We speculate on the missing ingredient 
that gives rise to that difference in Section ITVl 

A. Fractal Scale-Free Nets 

Consider the (2, 2)-flower, as a prototypical example 
of fractal hierarchical scale-free nets. In this net the dis- 
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tance between opposite hubs (or the diameter) scales as 
L n ~ 2™, and the mass scales as N n ~ 4™ ~ L\, cor- 
responding to a fractal dimension of d — 2. Suppose 
that each link is present with probability p. Denote the 
probability for two opposite hubs in generation n to be 
connected by P(p), then P', the analogous quantity in 
generation n + 1, is 



where 



pi = 2 p2 _ pi^ 



(8) 



Indeed, on iterating the flower to generation n + 1 the 
probability of contact between opposite hubs depends on 
the existence of either of two parallel paths, each consist- 
ing of two stringed copies of generation n. The proba- 
bility that one of the paths is connected is P 2 and ((SJ) 
follows: a naive addition P 2 + P 2 over-counts the event 
that all 4 generation-n units are connected (probability 
P 4 ). 

Starting with P = 1, for generation n = 0, one can 
then compute the probability of contact P for any other 
generation (Fig. . Eq. © has an unstable fixed point at 
Pc = (V5-l)/2 « 0.618 [H, where P'{p c ) = P{p c ) = p c , 
and two stable fixed points at P — and P — 1. If p > p c , 
the contact probability flows to P = 1 (as n — > oo) and 
the system is in the percolating phase. For p < p c , it 
flows to P — and there is no percolation. 
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FIG. 4: (Color online) Contact probability between opposite 
hubs in the (2, 2)-flower. Shown are the curves for generations 



Near the percolation phase transition the contact prob- 
ability obeys the finite-size scaling relation 



P = <E>( ApL 1 /" 



(9) 



where Ap — p — p c , and v is the critical exponent gov- 
erning the scaling of the correlation length, 

z,~\v-vc\- v . (10) 

We can obtain v by evaluating the derivative of Eq. © 
at p = p c - 



dP' 
dp 



A 



dP 
dp 



(11) 



A 



dP' 
~dP 



(12) 



Using ([9]) it then follows that 

A = {L n+l /L n fl v 



(13) 

6 — 2y5 and L n+ i/L n = 2, 



In our case A = 4(p c — p c 
yielding v = 1.63528 . . . 

Next, we address the probability that a site belongs 
to the infinite incipient cluster (or the giant component), 
Poo(p)i hi generation n. It obeys the scaling relation 



N- e ^(ApN e/p ). 



(14) 



The finite-size scaling exponent characterizes the size 
of the giant component at the transition point, p = p c : 



The scaling function has a non-analytic part 



^(x) ~ x@ , for small x, so that near the transition point 

Peo~(p-Pcf- (15) 




FIG. 5: Contact probabilities used in the derivation of 
Eq. |[TSJ). See text. 

Let A, P, C, and D denote the probabilities that a site 
is connected to exactly one, two, three, or four of the 
hubs, respectively (Fig. [5ji), then P QO =A + B + C + D. 
The analogous quantities in generation n + 1, are 



A' = SQ, 

B' = SPQ + TQ 2 , 

C = SP 2 Q + T(2PQ 2 ), 

D' = 5P 3 + T(3P 2 Q + P 3 ). 



(16) 



Here Q = 1 — P is the probability that opposite hubs (in 
generation n) are disconnected. S (T) denote the event 
that only one (two) of the hubs that the site reaches in 
generation n are also hubs of generation n + 1 (Fig. 0d). 
These are straightforwardly related to the A, B, C, D: 



S 



T 



-A + B 



:C, 



(17) 



-C + D. 
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As a useful check, one may verify that — A' + B' + 
C'+D' = S+T. From ([TB|) and $1T§ we obtain a recursion 
relation for S and T: 



S' 
V 



±Q(1 + P) 2 Q 2 (l + P) , 
1 -P 2 fl + P) P + P 2 - P 3 I ' ^ ' ^ ' 



The scaling of the giant component is dominated by A, 
the largest eigenvalue of the above matrix, evaluated at 

P = Pc, 



(N'/N)' 9 = A. 



(19) 



In our case A = (7- 2^5 + V73^32V5 )/4 and A'/A = 
4, yielding = 0.0503564 ... 

To obtain j3, we derive Eq. (TT4"]) with respect to p, 



0P* 



9p 



o 1-/9 

iV e — ^-tf(0) 



op 



dP 

Pa Op 



where we used the fact that Poo(0) ~ A™. Doing the same 
for P' and dividing the two relations, while using (jTTj) . 
we get 



(N'/N) 9 ^ 



AA, 



as n 



(20) 



Substituting for the values of A, A, N'/N, and 9, we find 
for the (2, 2)-flower j3 = 0.164694 ... 



Poo 0.6 




FIG. 6: (Color online) Poo for the (2, 2)-flower. Shown are 
curves for generations n = 6,7, 8, 10, and oo, obtained from 
Poc — S + T and iterating Eqs. ([T5J) and @. 



The scaling relations between critical exponents, fa- 
miliar from percolation in regular and fractal spaces, are 
obeyed as well. Indeed, we may rewrite Eq. (| 13jl as 

A = {N'/N) 1 '" 11 

since d is the fractal dimension of the hierarchical flower 
and N ~ L d . Using this, in conjunction with Eqs. fTJ 
and (|20p . we derive the scaling relation 



(21) 



o=4' 

vd 

The giant component, at criticality, scales as 
N g - L dg - N d o /d , 



where d g is its fractal dimension. Comparing this to 
N g ~ N x ~ 9 , on the one hand, and to Eq. (f2"Tj) , on the 
other hand, we get 



d g = d 



(22) 



which is a well-known scaling relation for percolation in 
regular space [2l|, [22j • 

The analysis carried out above for the (2, 2)-flower may 
be extended for other values of u > 1 and v. The recur- 
sion relation between opposite hubs in the general case 
is 



pi _ ptl i pv _ p 



u+v 



(23) 



As for the (2, 2)-flower, this has two stable fixed points 
at P = and 1, and an unstable fixed point whose loca- 
tion at P = p c may be computed numerically. One can 
then evaluate the correlation length exponent v, using 
Eqs. (|12[) . (|13[) and keeping in mind that L n+ \/ L n — * u 
in the thermodynamic limit. Results from such calcula- 
tions are shown in Fig. [7] There is general agreement 
with percolation in regular d-dimensional lattices, espe- 
cially for the particular case of v = u (note the analytical 
continuation to non- integer values of u and v) . 



B. Small-World, Assortative Nets 



In summary, percolation in fractal scale-free nets is 
very similar to percolation in fractal and regular lattices. 
As far as we can tell, the broad scale-free degree distribu- 
tion does not give rise to any kind of anomalous behavior, 
different from percolation in regular spaces. Note, for ex- 
ample, that the contact probability, P (Fig. [4]), and Poo 
(Fig. [6]) follow the same pattern as for percolation in 
regular spaces. Interestingly, the exponents v — 1.63528 
and /3 = 0.164694 that we find for the (2, 2)-flower (whose 
fractal dimension is d — 2) are not very different from the 
corresponding exponents in regular 2-dimensional space: 
v = 1.5076 and = 5/36 = 0.13888. 



We next turn to small world nets, and consider the 
(1, v)-flower as examples of this type of networks. The 
mass of the (1, v)-flower grows like N n ~ (1 + v) n , while 
the diameter increases only logarithmically, L n ~ In N n , 
making it a small world net of infinite dimensionality. 
As we shall shortly see, there is no percolation transi- 
tion, contradicting the finding for percolation in random 
scale- free nets d,|j,[5, 6]. This may be perhaps attributed 
to the fact that (1, i>)-flowers are quite strongly assorta- 
tive (highly connected nodes tend to be connected to one 
another), making them particularly resilient to random 
dilution. 
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FIG. 7: (Color online) Correlation length exponent v for frac- 
tal (it, u)-flowers (it > 1), plotted against their dimensionality 
d = ln(ii+ij) / In u (o). Shown are results for increasing u (bot- 
tom to top) and increasing v (left to right). The solid curve 
corresponds to v = u — a case that is close, numerically, to 
regular d-dimensional space (+). 

The recursion relation for the probability of contact 
between hubs in successive generations is now 

P' = P + P V - P v+1 , (24) 

which has an unstable hxed point at p = and a stable 
fixed point at p = 1. In other words, regardless of the 
dilution level, p, contact between hubs is guaranteed, in 
the thermodynamic limit of N — > oo. 

Let Ai (i = 1, 2, . . . , v + 1) denote the probability that 
a node is connected to exactly i adjacent hubs, in gener- 
ation n. The recursion relations for the analogous quan- 
tities in generation n + 1 are 

A[ = SQ , 

A' 2 = SPQ + TQ 2 , 

(25) 

A'i = SP^Q + T(i - 1)P'- 2 Q 2 , v ; 



A' v+1 = SP V + T{vP°- 1 Q + P v ) , 

where 

S = ( A x + A 2 + ■ ■ ■ + A v ) , 
v + 1 V / 

T = — |— (A 2 + 2A 3 + ■ ■ ■ + {v - l)A v + (v + 1)A V+1 ) , 
u + 1 V / 

(26) 

have the same meaning as for the (2, 2)-flower, in the 
previous section. Analysis of these equations reveals that 
there is a finite probability P^ for a site to belong to the 
giant component, at any dilution level p. For p small, 

p - - Ghml) ' (2?) 



and there is an essential singularity at p = 0. Practi- 
cally, though, it is impossible to tell whether P^ = or 
not, for p sufficiently small, and one could not rule out a 
percolation phase transition at some p c > based on a 
numerical study or on simulations alone (Fig. [8]). 




V 

FIG. 8: (Color online) P^, for the (1, 3)-flower. Shown are 
curves for generations n — 6,7, 8, 10, and oo, obtained from 
P = J^i A i and iterating Eqs. flM} - (J25J). Note that P x 
is indistinguishable from zero, in the scale of the plot, for p 
below about 0.2. 



C. Small- World, Disassortative Nets 

Having failed to find a percolation transition in the 
assortative (1, u)-flowers, we now turn to the (2, 2)-flower 
with a non-iterated link (Fig. [3]). The recursion relation 
for the probability of contact between hubs in successive 
generations is 

P' = \-{l-p){l- P 2 ) 2 . (28) 

Indeed, note that contact can be made through either of 
the two paths consisting of two stringed copies of gener- 
ation n (with probability P 2 , in either case) or through 
the non- iterated link (with probability p). The probabil- 
ity that none of these three parallel paths make contact 
is therefore (1 — p)(l — P 2 ) 2 , and P 1 follows. 

In the thermodynamic limit, P' — > P. It is easier to 
obtain P(p) implicitly, inverting (|28|) : 

1 

P=1 ~ (1 + P)(1-P 2 )' 

One can thus see that P(p) is double-valued, for p < 
5/32. A stability analysis reveals that only the lower 
branch is stable. For p > 5/32, the only available solution 
to d28j) is P' = P = 1. This solution is stable as well. 
Thus, P(p) has a discontinuity at p c — 5/32, where it 
jumps from P(p~) = 1/3 to P{p^) = 1, see Fig. |H1 

The recursion relations for the giant component are 
slightly more involved than in previous cases. We de- 
fine, as usual, A, B, . . . ,G as the probabilities that a 
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node reaches various hubs combinations in generation n 
(Fig. [TUkV We also denote by X the probability that, 
after embedding the n-th generation in generation n + 1, 
the node reaches only one of the hubs, connected to the 
non-iterated link. Similarly, Y is the probability that it 
reaches a single hub that is not connected to the non- 
iterated link, and Z the probability that it reaches both 
hubs (Fig. [TOb). We then have 

A' = XqQ , 

B' = YQ , 

C = {X + Y)qPQ + ZqQ 2 , 
D' = XpQ 3 , 

E' = YqP 2 Q + ZqPQ 2 , (29) 
F' = X{ P PQ 2 + Q[qP 2 + p(2PQ + P 2 )}} 

+ Y P PQ 2 + ZQ 2 {qP+p), 
G' = (X + Y)P[p(2PQ + P 2 ) + qP 2 } 

+ Z(l-qQ 2 ), 

where 

X = ^(A + C + E) 1 

Y = ^(A + C + E), ( 30 ) 
Z = D + F + G . 

Again, the fact that A' + B' H VG' = X + Y + Z 

confirms that the equations are consistent. 

The recursion equations simplify with the substitution 
S = X = Y and T = Z, leading to 

(S'\_f kQ(l + P) 2 \qQ 2 (l + P) \(S\ 
\T> ) {Pt+ P Q(1 + P) l-qQ 2 (l + P) ){T)- 

(31) 



The scaling of near the percolation threshold is dom- 
inated by the largest eigenvalue of the recursion matrix, 
A = (3 + V6)/6 (evaluated at p c ). From this we find the 
finite size exponent 9 = - In A/ In 4 = 0.0694207 . . . 

In order to compute the order parameter critical ex- 
ponent, (3, we must first replace the scaling relations 
suitable for fractals and regular spaces with relations for 
percolation in small world substrata, that are infinite- 
dimensional. Instead of ((9]) we write 

P = $(ApN x \ (32) 

where the original scaling argument L x l v has been 
changed to N x , obviating the question of diameter. Us- 
ing this and a similar argument to the one leading to (I2ip , 
we now derive 

x=~. (33) 

P 

Naively equating this relation to (|2"Tj) one obtains x = 
1/vd. This is, of course, meaningless, but makes some 
kind of sense: because the net is small world both its 
dimension d and the inverse of the correlation length ex- 
ponent \jv are infinite, but in such a way that their ratio 
yields a finite x. 

The exponent x is obtained in practice from the recur- 
sion relation for the contact probability, and using (|32[) : 
(N'/N) x = A = dP'/dP\ p=Pc . For the decorated (2,2)- 
flower we find x = 0, so that (3 — oo. Iterated curves of 
Poo (fc> r n ~ * °o) indeed show a transition at p c = 5/32 
with an infinite order parameter exponent (3 (Fig. Hip . 

IV. DISCUSSION 

We have studied the percolation phase transition in a 
class of hierarchical scale-free nets that can be built to 
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5/32 p 

FIG. 11: (Color online) Poo for the decorated (2, 2)-flower. 
Shown are curves for generations n = 6,7, 8, 9, 10, 11, and oo. 
Inset: Detail about p = p c , showing that Poo ~ (p — pc) 13 with 
/3 = co. 

display a large variety of structural properties and that 
can be analyzed exactly. When the scale-free nets are 
also fractal, that is, when the mass of the net increases 
as a power of its diameter, percolation is very similar 
to what is found in regular lattices. We do not see any 
specific signature that might be ascribed to the scale-free 
degree distribution. 

Percolation in small world hierarchical lattices is more 
exotic. In the (1, v)-flowers, we find that there is no per- 
colation phase transition: the system is always in the 
percolating phase, even as the bonds get diluted to con- 
centration p — > 0. This is in line with what is known 
for stochastic scale-free nets of degree exponent 7 < 3. 
However, for (1, v) -flowers the percolation phase transi- 
tion fails to appear even as v (and 7) increase without 
bound. To be sure, Poo(p) flattens more pronouncedly 
about the origin, reaching near zero probability at wider 
and wider regions of p, as v increases, but there is no tran- 
sition nevertheless — see Eq. (|27]). A possible cause for 



the exceptional resilience of (1, u)-flowers is their being 
strongly assortative, compared to stochastic and every- 
day life scale-free nets. 

We then studied percolation in the decorated (2, 2)- 
flower, a hierarchical scale- free net of degree exponent 
7 = 3 that is small world and disassortative. In this 
case there is a percolation phase transition at a finite 
Pc and the order parameter critical exponent charac- 
terizing the transition is [3 = 00. This agrees with the 
result for percolation in stochastic scale-free nets, that 
(3 = 1/(7—3) 0, since the degre exponent of the dec- 
orated (2, 2)-flower is 7 = 3. However, the finding in 
our case is generic: one can show that (3 = 00 for deco- 
rated flowers with other 1 < u < v values, independently 
of 7. The same is true for the (2, 2)-fiower decorated 
with two non-iterated links (connecting the two pairs of 
opposite hubs). The decorated (2,2)-fiower and similar 
constructs closely mimic everyday life stochastic scale- 
free nets (small world, disassortative, and high degree of 
clustering). Why is it then that they cannot reproduce 
a phase transition with finite [31 — Perhaps we are still 
missing out on some crucial structural property, common 
to everyday life stochastic networks. Another possibility 
is that it is a consequence of the hierarchical flowers be- 
ing finitely ramified (they can be disjointed by removing 
a finite number of nodes, regardless of the graphs' sizes). 
We do not know whether finite ramification is typical of 
everyday life networks. Finding out the answers to these 
questions will shed further light on the structure of the 
complex nets around us. 
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